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Abstract 

We identify the algorithm for constructing steady states of the n-species totally asymmetric 
simple exclusion process (TASEP) on L site periodic chain by Ferrari and Martin with a compo¬ 
sition of combinatorial R for the quantum affine algebra Uq^slif) in crystal base theory. Based 
on this connection and the factorized form of the R matrix derived recently from the tetrahedron 
equation, we establish a new matrix product formula for the steady state of the TASEP which is 
expressed in terms of corner transfer matrices of the g-oscillator valued five-vertex model at q = 0. 


1. Introduction 

In this paper and the next m we launch a new approach and integrable structure on the 
(n + l)-state totally asymmetric simple exclusion process on a one-dimensional (ID) periodic chain 
Tjl of L sites. It is called n-species TASEP or n- TASEP for short. A local spin at each site takes 
values in {0,1,..., n} and the neighboring pairs ( a , /3) with a > f) are exchanged to (/3, a) with a 
uniform probability. It is a model of non-equilibrium stochastic dynamics in physical, biological 
and many other systems including traffic flow, etc. See for example mm and reference therein. 

The first basic problem in the n- TASEP is to determine its steady state, which exhibits an 
intriguing non-uniform measure for n > 2 as in Example 13.11 It has been solved in m by 
introducing a companion system on an n-tuple of (0, l}-sequences called multiline process. Steady 
states of the multiline process possess a uniform measure and those for the n- TASEP are obtained 
as the image by a certain projection 7r from the former to the latter. The map 7r playing the key 
role is constructed via a combinatorial procedure on the n-tuple of {0, l}-sequences which we call 
the Ferrari-Martin algorithm. 

Our first main observation is that the Ferrari-Martin algorithm is most naturally formulated in 
terms of combinatorial R in crystal base theory, a theory of quantum group at q = o m- More 
specifically, each {0, l}-sequence in the multiline process is regarded as a crystal B l (12.11) of an 
antisymmetric tensor representation of the Drinfcld-Jimbo quantum affine algebra U q {slif). The 
“2-body operation” of the adjacent {0, l}-sequences described as arrival/service/departure [TT] is 
nothing but the combinatorial R acting on the tensor product B l ® B rn of crystals [21]. It is the 
bijection that arises from the quantum R matrix at q = 0 (crystallization) which still satisfies the 
Yang-Baxter equation. The “n-body operation” in the multiline process is identified, by using 
the Yang-Baxter equation, with the corner transfer matrix [3] of the size-n vertex model whose 
Boltzmann weights are crystallized to the combinatorial R. See (11.21) . We remark that such a 
size-n system equipped with the internal symmetry U q (slL ) corresponds to the cross-channel of 
the original n- TASEP in the sense that the role of the physical space Zz, and the internal space 
{0, ...,n} is interchanged. In particular the cyclic symmetry of the physical space Zz in the 
n-TASEP has been incorporated into the Dynkin diagram of U q (slL)- The integrability in the 
direct-channel, i.e. a relation to a U q (sl n +\) spin chain on Zz, has been demonstrated for a more 
general ASEP in [2] Sec.5]. 
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Our second result is the matrix product representation (14.411 of the steady state probability. It is 
expressed in terms of the operator Xi (14.611 which can be regarded as a corner transfer matrix of the 
five-vertex model (12.201) whose Boltzmann weights take values in a g-deformed oscillator algebra at 
q = 0 (12.161) . It is obtained by combining the crystal formulation of the Ferrari-Martin algorithm 
with the factorized form of the combinatorial R (12.221) . The latter is the q = 0 corollary of the 
factorized form of the quantum R matrix established recently [20j . The origin of the factorization 
itself goes further back to a 3D generalization of the Yang-Baxter equation known as the tetrahedron 
equation |24j . In fact our operator Xi , which is different from the earlier one m, possesses a far- 
reaching generalization in the framework of a 3D lattice model associated with the tetrahedron 
equation [T9] . 

The layout of the paper is as follows. In Section [2] we recall the background from U q (sli) and 
the combinatorial R necessary in this paper. In Section [3] we treat the n-TASEP and reformulate 
the Ferrari-Martin algorithm in terms of crystals and combinatorial R. In Section [4] a new matrix 
product formula for the steady state probability is derived by synthesizing the contents in Section 
[2] and Section [3] Section [5] contains a remark on a generalization of stochastic dynamics from the 
viewpoint of crystal base theory and a brief announcement on the so called hat relation. There are 
notable 3D structures connected to the tetrahedron equation behind several scenes in this paper. 
They will be fully demonstrated in our forthcoming paper M- 


2. Background facts from U q (sl L ) 


Let us recapitulate relevant results from the representation theory of quantum affine algebras, 
matrix product forms of the quantum R matrix and combinatorial R. 


2.1. U q (slL) and antisymmetric representation. The Drinfeld-Jimbo quantum affine algebra 
U q = U q (slL) is a Hopf algebra with generators ej,/i,fc ! ±1 (i £ Z l) [B; US] satisfying the Serre 
relation and 


kiki — k- ki — 1, \k t . kj] — 0, k'i e.j — q ■' ej kj . kifj — q ‘' 3 fjki, [e?;, fj] — Sij 


kj — k- 


q-q 


-i > 


e i e i ~(q + q )e i e j e i + ejej = 0, f, fj -{q + q )fifjfi + =0 (i - j = ±1 mod L). 


Here dij = 2 Sij —6i t j —i — Sij+i with dij = 1 if i = j modL and 5ij = 0 otherwise. All the indices 
are to be understood in Z l likewise. We adopt the coproduct A : U q —> U q ® U q of the form 

A et = et <g> 1 + ki <8> e i; A f z = 1 ® /) + /i <8> kf 1 , Afc^ 1 = kf 1 <g> kf 1 . 


Our A here is opposite of [2D:. We assume that q is a generic complex parameter in Section liPl 12.31 
We are concerned with the degree l antisymmetric tensor representation <p x : U q {sli,) —>• EndH* 
with 1 < l < L. The representation space V 1 is given by 

y*=0C|b>, B l = {h= (6r,..., 6 l) e {0 ,\} L \ |b| = l}, (2.1) 

b GB l 

where |b| = b\ + ■ ■ ■ + • Thus dim W = (^). The labeling set B l of the basis of V 1 is called 
the crystal of V 1 . Note that the dependence on L is not indicated in the notation B l . To save 
the space, elements of B l will often be written as words on {0,1}, .e.g. (1, 0,1,1, 0) £ B 3 is just 
denoted by 10110, etcQ. 

Generators act on V 1 as ((j> x (g) for a generator g is denoted by g for simplicity) 

ei|b) = x 5i, ° |b + e; - e i+ i), fi\b) = a; _<5i ' 0 |b - e* + e i+ i), ^ ±:L |b) = q M b i~ b i+i) |b), 


where ei = (0,..., 0,1, 0,..., 0) £ Z L and the right hand sides are to be understood as 0 unless 
the arrays in the ket stay within {0,1} L . The representation <f> x is irreducible and x is called a 
spectral parameter. 




more standard notation is a one-column semi-standard tableau, the one filled with 1,3,4 for this example. 
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2.2. Quantum R matrix. Let A X:V = {<j) x <g> <j > y )A be the tensor product representations on 
V 1 8 V m , where 1 < l,m < L are arbitrary. It is known that they are irreducible if z := x/y 
is generic. Moreover there is a linear map lR(;r) = J{ l ’ m (z) : V 1 8> V m —> V m 8> V 1 called a 
quantum R matrix which is uniquely characterized by the intertwining relation with the quantum 
group A Xj y(g)tR(z) = Jl(z)A yiX (g) for g £ U q (sIl) up to an overall scalar. Let us write its matrix 
elements on the basis |i) <g> |j) 6 b 1 ® V m as 

*(*)( |i>®tf>)=£*(*)?j b |b>®|a), (2-2) 

a,b 


where (|i|, |j|) = (Z,m) and the sum extends over a, b such that (|a|, |b|) = (l,m). The matrix 
elements have the support property reflecting the weight conservation: 

lR(z)fj b = 0 unless a + b = i+j (e {0,1, 2}^). (2.3) 

The R matrix satisfies the Yang-Baxter equation j3] 

(3l'- m (s) 8> 1)(1 (g) ^{zz'))^' 1 ^') ® 1) = (1 ® ft fe ’V))(^ fc ’ m (-^') ® 1)(1 ® % l ’ m {z)) (2.4) 

for any 1 < k. I, m < L. We normalize the Jl l,m (z ) so that 

min(Z,m) — 1 

■^(~)(l e <j) 8> |e< TO >) = p(z)|e< m ) 8> |e<j), g(z) = (1 - (-1 ) l+rn q l+m ~ 2l z), (2.5) 

L) + 

where (m)+ = max(m, 0) and |e<j) = |ei + • • • + e;) £ V 1 . Then all the elements lR(*)fj b are 
polynomials in q and z. 


Example 2.1. The IR 1,1 ^) is the well known R matrix for the vector representation of U q {slL,)- 
The L = 2 case corresponds to the six-vertex model [3]- The nonzero elements read 

a(*)2£ = i - 9 2 *, = 9(1 - z), = (1 - ?V (i< ^ 

where 1 < i ^ j < L and (9(true) = 1 , 0 (false) = 0 . 


Example 2.2. Nonzero elements of fR 1,2 ^) for U q (sl 3 ) are as follows. 

rpioo.no _ rpOio.no _ rpioo.ioi __ rpooi.ioi _ rpoio.on __ rpooi.on _ -1 , „ 3 „ 

*100,110 — *010,110 — *100,101 — *001,101 — -*010,011 — -*001.011 — 1 + 9 *> 

rpOOl.llO _ rp 010.101 _ rp 100,011 __ /. \ rpOOl.llO _ rp 010.101 _ rplOO.Oll _ _ /-t _ 2 \ 

*001,110 — *010,101 — *100,011 — 91 1 ' 9*ji *010,101 — *100,011 — **001,110 — 9 U 9 J z > 

rpOOl .110 _ rpOlO .101 _ rplOO.Oll _ /I 2 \_ 

*100,011 _ ■‘'*001,110 — **010,101 ~ v 1 9 ) z - 


Example 2.3. Nonzero elements of CR 2,1 (z) for U q (sl 3 ) are as follows. 

rp 110.100 _ rp 110,010 _ rp 101,100 __ rpl01,001 _ rpOll.010 __ rpOll.OOl _ 1 , -3 _ 

*110,100 ~ *110,010 — *101,100 — *101,001 — *011,010 — *011,001 — 1 + 9 *> 

rpOll, 100 _ rp 101,010 _ rp 110,001 __ /, , \ rpOll.lOO _ rpl01,010 _ rpllO.OOl 

* 011,100 — * 101,010 — * 110,001 — 91 1 + 9*J> * 101,010 — * 011,001 — ** 011,100 

rpOll,100 _ ,rpl01,010 _ rp 110,001 _ fl _ -2\„ 

*110,001 _ "*011,100 — **101,010 — V 1 9 I*. 


-9(1 -9 2 )*, 


See for example [ 6 ] for more information on the quantum R matrix. 


2.3. Matrix product representation of 3t(z). In a recent work m based on the tetrahedron 
equation, a matrix product representation of the quantum R matrix 3i l ’ m (z) was constructed. Let 
us quote the result in a form adapted to the present convention. 

Consider the g-oscillator algebra A q generated by a + , a - , k satisfying the relations 

ka ± = -q ± 1 a ± k, a+a" = 1 - k 2 , a”a+ = l- g 2 k 2 . (2.6) 

They may be identified with the operators acting on the Fock space F = © m>0 C|m)0 as 

a + |m) = \m + 1), a~ \m) = (1 — q 2rn )\m — 1), k|m) = (— q) m \m). (2.7) 


2 This ket | m) £ F with m € Z>o should not be confused with |b) £ V 1 with b (= B 1 . 



4 


ATSUO KUNIBA, SHOUYA MARUYAMA, AND MASATO OKADO 


We arrange them in the 3D L operator £ = (£©) [4j where all the indices belong to {0,1}. They 
are zero except the following: 


n 0,0 n 1,1 -i /"»U,1 

^ 0,0 — ^ 1,1 — ^ 1,0 — a 


'1,1 


0,1 


r 1 ’ 0 - « 
^0,1 — a 


L 


0,1 

0,1 


= <l k - 


The £ may be viewed as defining a six-vertex model whose Boltzmann weights belong to A 
q-oscillator valued six-vertex model: 


( 2 . 8 ) 

g , i.e., 


n a,b _ 

L i,0 ~ 



0 1 


1 


0 1 0 



1 


1 1—► 0 0 


1 0 


0 1 


1 


0 1 
1 1 


Oil 
a + a~ k 


0 

qk 


(2.9) 


The relations in (USD involving a ± a =F are quantization of the so called free fermion condition 
[3J eq.(10.16.4)]| aJ7=aJs= o. One may introduce the third arrow perpendicular to the two arrows 
for each vertex encoding the action of the corresponding g-oscillator generators. It amounts to 
regarding £ as a vertex of the 3D square lattice, hence the name 3D L operator. The resulting 3D 
classical spin model satisfies the tetrahedron equation [201 eq.(2.13)] involving the 3D R operator 
32 obtained by replacing q by — q in [201 eq.(2.2)]. See also [3]. Such a 3D structure is a source of 
the factorization of the R matrix in Theorem 12.41 and ultimately the matrix product form of the 
steady state probability that we will derive in Section [I] 

Let h be the operator on F defined by h|m) = m|m). Denote by F* = © m> 0 C(m| the dual of 
F with the pairing specified by (m|m / ) = ( q 2 ) m 5m,m' with ( q 2 ) m = TljliCl — g 2 © Accordingly we 
set Tr(z h X) = J2 m >o for X e Ar Set e(z) = g- (i " m) +( 1 - (~l) l + m z)g(z) where 

p(z) is defined in ([2.51) . The following result is a special case ei = ■ • • = £l = 1 of [201 Th.4.1], 


Theorem 2.4. Elements of the R matrix 12. £1) with the normalization JO) are expressed in the 
matrix product form 


a,b 




p(z)Tr 


(z h L 


ai,bi 

*liJi 


r a L ,b L \ 

^iL,jh ) ’ 


where a = (ai,..., a^) and similarly for b, i, j. 


( 2 . 10 ) 


For instance we have 


rnOl 1,100 
^ 011,100 


e( 2 )Tr(z h £° ; i£[ ; °£[;°) 


g 2 p(z)Tr(z h k 3 ) 


g 2 e(z) 

1 + q 3 Z 


in agreement with Example 12.31 due to (l,m) = (2,1) and g(z) = g _1 (l + qz){ 1 + g 3 ^). 

We remark that the operators e, 6, A used for the n-ASEP [22] can be realized in A q 1/2 (12.61) via 
ei-)a + ,^a“,dH'k 2 . 


2.4. Combinatorial R. The product set of crystals B l and B m is again called a crystal and 
denoted by B l ® B m . As mentioned before, all the elements of the quantum R matrix = 3i l ' m : 
V 1 V m —>• V m 8) V 1 are polynomials in q. Therefore it makes sense to set q = 0. We thus 
define R = R l ’ m = 3l l ’ m {z = 1)| 9= <11. It turns out that R becomes a beautiful bijection among the 
crystals R l ’ m : B l 8 1 B m —> B m 8> B l , which is called a combinatorial R. 

Example 2.5. Example 12.21 with (12.21) yields the combinatorial R 1 ’ 2 : B 1 8 > B 2 —>• B 2 81 B 1 as 
100 8) 110 i-a 110 8> 100, 100 8) 101 ^ 101 8) 100, 100 8) 011 1-4 110 8) 001, 

010® 110 i-a 110 8>010, 010 ® 101 ^ 011® 100, 010 ® Oil ^ 011® 010, 

001® 110 I-A 101 ® 010, 001® 101 HA 101® 001, 001 ® Oil I-A oil® 001. 


3 Although z will not be concerned in this paper it is an important ingredient of the combinatorial R in the 
literature. 
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Example 2.6. Example 12.31 with (12.211 yields the combinatorial R 2 ' 1 : B 2 ® B 1 M B 1 ® B 2 as 

110® 100 m 100® 110, 110® 010 m 010 0 110, 110® 001 m 100® Oil, 

101 0 100 m-100 0 101, 101 0 010 m 001 0 110, 101 0 001m 001 0 101, 

011 0 100 m 010® 101, 011 0 010 m 010 0 011, 0110 001m 001 0 011. 

The maps R 21 and R 1,2 are both bisections inverse to each other. They are nontrivial in the 
sense that R( i 0 j) yf j 0 i in general. 

Back to the general case, let R = R l,m : B l ® B m M B m ® B l be the combinatorial R. Write 
its matrix element like (El: 

R( j 0 j) = Rf ; b b ® a. (2-11) 

a,b 

Then for any i 0 j, the matrix elements f??j b are all 0 except exactly one pair b ® a G B m ® B l 
determined from i and j. An algorithm for finding it, namely the image b ® a = R(i 0 j), was 
obtained in [23 Rule 3.100. We call it NY(Nakayashiki-Yamada)-rule. 

NY-rule for l < m: 

(i) Given i ® j G B l 0 B m , draw a tableau for i below and the one for j above. For example, the 
tableaux for i ® j = 1100100100® 0010111110 G B 4 ® B 6 look as (i) below (1 and 0 are denoted 
by • and empty box): 


(i) (ii) (ih) 



(ii-1) Pick any dot (call it d) in i and connect it to a dot (call it d') in j. The partner d' must be 
the rightmost one among the dots located weakly leftQ of d. If there is no such dot, return to the 
right (periodic boundary condition) and d! is the rightmost one. 

(ii-2) Repeat (ii-1) for all the remaining unconnected dots in i. 

(iii) The image b ® a is obtained by shifting the (m — l) unconnected dots in j to the bottom. In 
the above example, b ® a = 1110110100 0 0000101110 G B 6 ® B A . 

Remark 2.7. In the step (ii), the order of picking unconnected dots in i is not unique, but the 
final result is independent of it due to [23 Prop. 3.17]. 

In the above diagram (ii), the dots in the bottom tableau are picked for example in the order 
2,1,3,4 from the left. The connecting lines are called R-line. We remark that NY-rule for R : 
B l 0 B m M B m ® B l for l < m implies 

R( i ® j) ^ b ® a => i < b and a < j, (2-12) 

where x < y is defined by y — x G (Z> 0 ) L for x, y G Z L . 

The NY-rule for R l ’ m with l > m is a dual of the above. One starts from a dot d in j G B m 
(upper tableau) and seeks a partner d' in i G B l (lower tableau). The d' should be the leftmost dot 
among those located weakly right of d in a periodic sense. By construction R rn ' l R Lj n = 
holds for any l and m. Note that R l ’ m is trivial at l = m in the sense that R l ' l ( i ® j) = j ® i in 
agreement with Example 12.11 with (q, z) = (0,1). 


^ Actually we explain the algorithm for its inverse since the assumptions l < m here and k > l in mi corresponds 
to the opposite situation. 

^ “Weakly left” means exactly above or strictly left. “Weakly right” means similar. 
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It is customary to depict the relation R( i (g> j) = b <g> a as 



(2.13) 


This formally looks same as m ■ Note however the arrows there carry 0,1 while those here do 
the elements from crystals B l ,B m which are L-tuples of 0,1. 

As the q = 0, z = z' = 1 corollary of (12.411 the combinatorial R also satisfies the Yang-Baxter 
equation: 

(R ® 1)(1 ® R)(R ® 1) = (1 ®R)(R® 1)(1 ® R). (2.14) 


For example the action of the both sides on 0100 ® 0011 ® 1101 G B 1 <8> B 2 <g> B 3 leads to 


0111 

1100 

0001 

0111 

1100 

0001 

X 

t 

t 

X 

0110 

1101 

0001 

0111 

1000 

0101 

t 

X = 

= X 

t 

0110 

0001 

1101 

0100 

1011 

0101 

X 

t 

t 

X 

0100 

0011 

1101 

0100 

0011 

1101 


Here = means that starting from the same bottom line one ends up with the same top line despite 
the different order of applications of the combinatorial R’s. Combinatorial l?’s constitute the 
most decent and systematic examples of set theoretical solutions of the Yang-Baxter equation [9] 
connected to the crystal base of quantum groups, which have numerous applications m\mm 

mmm- 


2.5. Matrix product representation of combinatorial R. Setting q = 0 in (12.101) leads to a 
matrix product representation of the combinatorial R. Let us write it out in terms of the 3D L 
operator and g-oscillator at q = 0. We will be exclusively concerned with R l,rn with l < m. 

First we define the q = 0-oscillator Aq to be the algebra generated by a + ,a~,k satisfying the 
relation^ 

k 2 = k, ka + = 0, a“k = 0, a a + = 1, a + a - = 1 — k. (2.16) 

They may be identified with the operators acting on the Fock space F = ©„ l>0 C|m) as 

a + \m) = \m + l), a - \m) = (1 - c> m ,o)|ro - 1), k|m) = <5 m ,o|m). (2-17) 

As a C-vector space, the Ao has a Poincare-Birkhoff-Witt type basis: 

1, (a+r, (a -y, (a+) s k(a-)‘, (2.18) 


where r > 1 and s,t > 0. We introduce the 3D L operator at q = 0 denoted by L = (L^j 1 ) with 
indices from (0,1}. Explicitly they are zero except the following: 


.0,0 _ r i,i _ . 

^ 0,0 — ^ 1,1 — 


7-0,1 _ 1 7 - U, J- _ 

J-jr\ 1 — K, -1 n — 3- 


^ 0,1 


0,1 


1,0 


7 1 ’ 0 - a- 
^0,1 — a 


(2.19) 


T a.b 

= * 


■a 0- 


■0 1- 


1 1- 


■0 0 - 


■1 0- 


1 

a 


1 

k 


( 2 . 20 ) 


° Although the notations aA k, F. etc are retained, they are to be distinguished from those for A q . 
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Note that compared with £ (12.81) . the “vertex” (the rightmost one in (12.91) i is missing due to 
q = 0. As the result L is regarded as a q = 0-oscillator valued five-vertex model whose Boltzmann 
weights belong to Aq. 

Let F* = © m>0 C(m| be the dual of F defined by (m|m'} = Let Ag n C Aq be the vector 

subspace spanned by (12.181) except 1. Then Tr(X) := X)m>o( m l^l m ) fi n i te f° r an y X G Aq 11 . 
From Theorem 12.101 and these definitions we have 


Corollary 2.8. Elements of the combinatorial R = R 1, 
are expressed in the matrix product form 


B l <gj B n 


Rtf = Tr (L 


ai,bi 

*l ,3 l 


. j^a L ,b L 
Wi JL 


)• 


B m <g> B l \2.11\) with l < m 


( 2 . 21 ) 


For example -Roiopoi = Tr(Lj’°I/5’jLg’}) = Tr(a~a+k) = 1 in agreement with Example 12.51 
From (12.191) the number of k’s contained in the above product is given by |b| — |i| = |j| — |a| = 
m — l > 0. Therefore Lf’f ''' Iff f' G Ag 11 is guaranteed and the trace is convergent. The absence 
of the vertex if q in (12.191) reflects the fact that there is no dot going up in the two row diagram in 
the NY-rulc of R l,m with l < m. In the product (12.211) . one can interpret that a + creates (emits) 
an ff-line from a dot in the lower tableau and a” annihilates (absorbs) an jT-line into a dot in the 
upper tableau. The state |to) G F corresponds to the “segment” where there are m H- lines. 

The L operator Iff (12.191) and the matrix product form of the combinatorial R (|2.21l) may be 
depicted in the 3D picture as follows: 


T a ->b _ 

Hi - 




( 2 . 22 ) 


Here the black arrows assigned with the indices carry 0 or 1. The blue ones carry the Fock space 
F over which the trace is taken. 


3. 71-species TASEP and its steady state 

3.1. n- TASEP. Consider the periodic ID lattice with L sites which will be denoted by Each 
site * G Zi is assigned with a physical variable (local state) cr,; G {0,1,... ,n}. It is interpreted 
as the species of the particle occupying it or 0 indicating the absence of particles. We assume 
1 < n < L throughout. Consider a stochastic model on such that neighboring pairs of local 
states (cr, ct') = (cr^, ) are interchanged as a a' —>• cr' a if cr > a' with the uniform transition 

rate. The whole space of states is given by 

(C n+i )0 L^ 0 C\a 1: ...,a L ), (3.1) 

(tri,...,erz,)e{0,...,n} i 

where we suppose that the ket | • • •) here can safely be distinguished from |b) G V 1 for b G B l in 
Section O and also from | m) G F for m G Z>o in m by the context. Let P(ai,..., <Jl', t) be 
the probability of finding the configuration (or,..., ol) at time t, and set 

I Pit)} = P(a 1 ,...,a L -,t)\o 1 ,...,a L ). (3.2) 

(a 1 ,...,<7 L )e{0,...,n} L 

By n-species TASEP, or n- TASEP for short, we mean the stochastic system governed by the 
continuous-time master equation 


j t \P(t))=H\P(t)), 


( 3 . 3 ) 
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where the Markov matrix (also called “Hamiltonian” by abuse of terminology despite it is not 
Hermitian in general) has the form 


E 


H — ^ ^- 2 , 2+1 5 


h\a, a') = 


\a',(j) -\a,a') (cr > o-')> 

0 (cr<a’). 


(3.4) 


Here hi^ + \ is the local Hamiltonian that acts as h on the ith and the (i + l)th components and 
as the identity elsewhere. As H preserves the particle content, it acts on each sector labeled with 
the multiplicity m = (mo,..., m n ) G (Z>o) n+1 of the particle^]: 


v(™)= y 

<rG5(n 


C|cr), S(m) = {a= (tJi,..., <jl) G {0, ...,nY 

) 


L 

E* 

3 =1 


'k,ui = mfc,Vfe}. (3.5) 


It has the dimension dim V (m) = llln! L, m r - The space of states is decomposed as (C n+1 )® L = 

©m=(m 0 m n ) ^( m )) where the sum ranges over nn G Z> 0 such that mo H-bm„ = L. A sector 

H(mo,... ,m n ) such that m, > 1 for all 0 < i < n is called basic. Non-basic sectors are equivalent 
to a basic sector for n'-TASEP with some n' < n by a suitable relabeling of species. Thus we shall 
exclusively deal with basic sectors in this paper (hence n < L as mentioned before). For various 
symmetry and spectral property of H, see [2]. In each sector V(m) there is a unique state |P(m)) 
up to a normalization, called the steady state , satisfying H\P(m)) = 0. 


Example 3.1. Up to an overall normalization, the steady states in V(l, 1,1), V(2,1,1) and 
V(l, 2,1) for 2-TASEP are given by 

|P(1,1,1)) = 2|012) + |021) + 1102) + 21120) + 2|201) + |210), 

|P(2,1,1)) = 3|0012) + |0021) + 2|0102) + 3|0120) + 2|0201) + |0210) 

+ 11002) + 211020) + 3| 1200) + 3|2001) + 2|2010) + |2100), 

|P(1, 2,1)) = 2|0112) + |0121) + |0211) + 11012) + |1021) + |1102) 

+ 211120 ) + 211201 ) + 11210 ) + 212011 ) + | 2101 ) + | 2110 ). 

Similarly the steady state in V(l, 1,1,1) for 3-TASEP is given by 

|P(1,1,1,1)) = 9|0123) + 3|0132) + 3|0213) + 3|0231) + 5|0312) + |0321) + 3|1023) + |1032) 

+ 5|1203) + 9| 1230) + 3|1302) + 3|1320) + 3|2013) + 5|2031) + |2103) + 3|2130) 

+ 9|2301) + 3|2310) + 9|3012) + 3|3021) + 3|3102) + 5|3120) + 3|3201) + |3210). 

One can observe the symmetry under the cyclic shifts which holds in general. 

A combinatorial algorithm to construct the steady state |P(m)) G V (m) for general basic sectors 
U(m) of the n-TASEP on Z l was obtained by Ferrari-Martin [lTj. In the rest of this section we 
show that it is most naturally formulated in terms of the combinatorial R in the previous section. 


3.2. Multiline process. To a basic sector U(m) of the n-TASEP with the multiplicity m = 
(mo,..., m n ), we associate another system called multiline process. Define the integers Iq, ..., l n +i 
by 

li = m n _ i+ 1 H-1- m n _i + m n (0<*<n + l). (3.6) 

They satisfy 0 = lo < li < ■ ■ ■ < l n < ln+i = L due to m, > 1 for all i. Given L, the data (mi) 
and (li) can be transformed to each other uniquely. Introduce the sets 

B( m) = B h ® ■ • • (g) B ln = (bi ® • ■ ■ ® b„ | bj G B lj }, 

B + ( m) = {bi ® ■ ■ • ® b n G B( m) | bi < • • • < b„} c B( m), 


^ V (m) here should not be confused with the antisymmetric tensor representation V m from GO- 


(3.7) 

(3.8) 
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where < is defined after (12.1211 . The following bijection is elementary and will be utilized later: 


V : 

S{ m) -A B+{ m); 

cr = 

(cti, ... ,a L ) ha y>i(er) 0 ■ ■ • 0 ^„(cr), 

(3.9) 

Vi 

(cr) = (9(a i >n+l 


• i 9(aj, >n+l— j)) £ B lj , 

(3.10) 

v~ 

^(bi 0 • ■ ■ 0 b„) = 

= bi + ■ 

1 • • + b n , 

(3.11) 

where 9 is defined i 

n Example 12.11 Applying 

• o -1 to (13.91) and usine (13.111) we see that 



o 

' = Vi( 

Cr) H-h 

(3-12) 

holds identically. If 

■ (13.101) is riven as 

Vi ( <T ) 

= ( 7 j.i,... ,7 j,l), it is also easy to see that 



cr k = i => 

(7i, k, ■ ■ 


(3.13) 


for each k £ . The sum in (13.111) as elements in Z 1, is not an operation usually performed in the 

crystal base theory. 

Example 3.2. Let n = 3 and take <r = (3, 0, 1, 2,3, 0,1) £ 5(m) with m = (2, 2, 1, 2). The image 
ifi(cr) = (fi(cr) <8> </? 2 (cr) 0 <^3 (o') £ £>+(m) C B 2 ® B 3 ® B 5 is given by 


cr 

3012301 

(Tk 

V^) 

1011101 

73,fc 

V2(cr) 

1001100 

72, k 

wW 

1000100 

7i,fc 


One can check tpi(<r) < ^{(t) < <p 3 (cr) (13.81) . cr = <^i(cr) + Vi{(t) + V:i( (T ) <13.1211 and (13.1311 . 

The multiline process eh is a stochastic process on the set £>(m). To describe its dynamics we 
first introduce a deterministic map 

T : B( m) <g> Z L -A Z L 0 B( m); bi 0 • • • 0 b„ <8> k ha k' 0 b( 0 • • • 0 b'„. (3-14) 

It is defined by sending Z^ from the right to the left through B l 1 <g> • • • (g) by successively 
applying the following pairwise exchange rule: 

B l < 8 > Z L — > Z L 0 B l ; (xi, ..., Xl) 0 A: ha k! 0 (x ' l5 ..., x^), 

min(xfc, Xfc+i) i = k, (3 15) 

max(x/c, Xfe+i) i = k+ 1 , 

Xj otherwise. 

The map T in (13.1411 induces the “time evolutions” on B(m) by setting 

T k : B(m) -A B(m); T fc (bi 0 • • • 0 b„) = bi 0 • • • 0 b^ (ft £ Z L ), (3.16) 

where the right hand side is the one appearing in (13.1411 . 

Example 3.3. Take s = bi 0 b 20 b 3 = 00001000010100001011 £ B 1 ®B 2 ®B 3 . ThenT 2 (s) = s 
and the other cases are given by 


k' = k + x k — 1, x' = < 


002031 002031 000231 002031 002013 102030 

t>3 

b 2 

bi 


Here the same tableaux as for the NY-rule are used. As an example, T$( s) has been obtained by 
applying (|3.1511 successively as 

s 0 3 = 000010 0 001010 (8> 001011 0 3 ha 000010 0 001010 <S> 3 0 000111 

HA 000010 ( 8 ) 3 0 000110 0 000111 ha 2 0 000010 0 000110 0 000111 = 2 0 T 3 (s). 


Ti(a) 


T 3 ( s) 


T 4 (s) 


Tb( s) 


T 6 ( s) 
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The sequences above the tableaux will be explained after (13.1911 . 

Now we can define a stochastic dynamics on B( m) by declaring that each element of it undergoes 
the L kinds of time evolutions Tj,.... TV, with an equal probability. The resulting system is called 
multiline process m- We have not found a crystal theoretic interpretation of Tfc. See Section PT~ 2 l 
for a further remark. 

3.3. Projection from multiline process to n-TASEP. Fix L,n and m = (m^) and (If) as 
explained around (13.61) . Recall that the sector V (m) of the n-TASEP is defined in (13.51) and the 
set of states of multiline process is £>(m) = B h <g) ■ • • <g> B ln in (13.71) . We introduce a map 

7 Tj : B( m) —> B lj ] s = bi <g) • • ■ <g> b n -jTj (s) (1 < j < n) (3-17) 

by a composition of the combinatorial R. It is most easily grasped by the diagram 

bj---^ 7 Tj(s), 

bj+i bj +2 b n (3.18) 

which is a composition of (12.131) . Namely one lets the component B lj penetrate through its 
right neighbor B lj + 1 <g) ■ ■ ■ <g> B ln by successively applying the combinatorial R = R l i' lr for r = 
j + 1, • • • , n. Note that nj (bi <g) • • ■ (g) b„) actually depends only on bj <g> ■ ■ ■ (g> b„ and in particular 

7r n (bi <g> • ■ ■ ® b„) = b n . _ 

From 1 1 < • • • < l n (see around (USD) and (12.121) . we have an increasing sequence 7Ti(s) < 
tt 2 (s) < • • • < 7r n (s) of elements of the crystals B l1 , B l2 ,..., B ln . In fact one can check 

7Ti(s) (g) • • • (g> 7r n (s) S B+(m) (3.19) 

for any s E B{ m). Sending this by ip^ 1 and applying (13.111) . we are able to define a surjective map 
7t from the states B( m) (13.71) of the multiline process to the states 5(m) (13.51) in the n-TASEP by 

7T : B(m) ->■ S'(m); 7r(s) = 7Ti(s)-f - f7r„(s). (3.20) 

The sequences above the tableaux in Example 13.31 show the images in 5(3,1,1,1) under tt. By 
comparing the NY-rule and the diagram (13.181) with Fig.2 and Fig 3 in [IT] and noting Remark 
12.71 it can be seen that the map 7r coincides with the Ferrari-Martin algorithm V^ : X —> V in 
eh upon conventional adjustment. This is our first main observation in this paper. 

As a consequence of this identification, Theorem 4.1 and Theorem 2.2 in EH lead to Proposition 
13.41 and Proposition 13.51 given below. 

Proposition 3.4. Let Ti be the map on S( m) that changes the ith and ( i+l)th components (cq, Ui+i) 
of cr E 5(m) into (min(<7j, <7i+i), max(cq, <7j+i)). Then the following diagram is commutative: 


B( m) ———> B( m) 



£(m) ———> 5(m). 

One can easily check the commutativity in Example 13.31 See Section [5721 for a crystal theoretical 
generalization of Ti and Ti. 

Proposition 3.5. Up to a normalization, the steady state in the sector V (m) is given by 

l-P( m )) = W 6 *)) = P ( cr )l <T )> 

sG^(m) <r=(tri ,...,cr^)G5(m) 

P(cr) = #{s E B( m) | 7r(s) = cr}. 


(3.21) 
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The last condition is equivalent to ip 1 (711 (s) ® ■ ■ ■ <g> 7 r„(s)) = er due to (Id. Ill) and (Id. 201) . By 
taking ip further and using (ld.Q|l we find 

P(er) = #{s e B{ m) | 7r fc (s) = tp k (<r) (1 < k < n)}. (3.22) 

In the next section we invoke the results in Section [53] to derive a new matrix product formula for 
the (unnormalized) steady state probability P(cr) of the configuration er. 


4. Matrix product formula 

4.1. Diagram for 7Tj: Combinatorial corner transfer matrix. There is a single diagram that 
represents all of 7 Ti(s),..., 7r„(s) in (13.181) simultaneously. We illustrate it along n = 3 case. Given 
s = bi b 2 ® b 3 € B( m) = B l1 Cg> B l2 (g) B ls , the diagrams (13.181) for 7Ti and 712 are 


bi 


b 2 


7Tl(s) 

b 3 


b 2 


7T 2 (S) 


b 3 


Here each vertex stands for the combinatorial R. See (12.131) . This defining diagram of 7 Ti(s) can 
be deformed by means of the Yang-Baxter equation (12.141) : 


bi 



7Tl(s) 


bi- 

>< 


7Tl(s) 


b 2 b 3 b 2 b 3 

Inclining the vertices appropriately in the last diagram we find that 7 Ti(s), 7r 2 (s), 7 t 3 (s)(= 
all be included in the left diagram below: 


n = 3 


n = 4 


(4.1) 
b 3 ) can 


bi 

b 2 

b 3 




t 



7T 2 (S 

S ) 


7T 3 ( 


bi 

b 2 

b 3 

b 4 


-t 7Ti(s) 


7T 2 (S) 


7T 3 (s) 


7T 4 (s) 


(4.2) 


The result for n = 4 case is also given together for which s should be understood as s = bi (g> 
• • • <g> b 4 e B( m) = B l1 ® • • • <g> B' 4 . General case is clear and similarly obtained by transforming 
the defining diagrams (13.181) of 7Tj(s) for j = 1 , ,n into composable forms with the aid of the 
Yang-Baxter equation as in (ED- 

In (14.21) the arrow starting from b^ carries the crystal B li and each segment corresponds to an 
element of it. The only nontrivial events are the combinatorial R's indicated by the vertices. The 
90° left turns along the arrows do not change the elements. It is consistent with the bottom line 
since 7r n (s) = b„ as noted after (13.181) . The inclusion of such turns makes the diagram symmetric in 
that there are n incoming arrows horizontally and n outgoing ones vertically. A further significance 
will become manifest in m- Since the combinatorial R is a deterministic map, every part of the 
diagram is fixed uniquely upon choosing the input s = b 4 <8> • ■ ■ ® b„ G B{ m) from the left. 

In this way we are led to the corner transfer matrix [3l Chap. 13] of the vertex model associated 
with the antisymmetric tensor representations of U q {sli,)- It is at the combinatorial point q = 0 
where every vertex is crystallized to the combinatorial R. Note that the physical space 'Ll and 
the internal degrees of freedom {0,1,..., n} in the original n- TASEP have been interchanged; we 
have a corner transfer matrix on the system of linear size n whose local interaction is encoded in 
the quantum group U q (slL) at q = 0. We stress that such a cross-channel of the problem has been 
captured precisely by reformulating the Ferrari-Martin algorithm in terms of combinatorial R and 
thereby enabling a systematic use of the Yang-Baxter equation. 
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4.2. Factorization into five-vertex model. By the graphical representation (14.211 of 7Tfc(s), the 
stationary probability (13.2211 for the configuration er = (ay,... , cjl ) is expressed, for instance for 
n = 3, as 


P(<T) 


E 

bi®b2®b 3 eB(m) 


bi ■ 

b 2 

b 3 — 1 


J 




<P2 (O’) 


<^3(cr) 


(4.3) 


Here each ipj{cr) signifies the boundary condition 7 Tj(bi (g> b 2 <8> b 3 ) = tpj(<r) that the element 
7Tj (bi <8> b 2 <S> b 3 ) € B lj making the left turn there should coincide with ifij(cr) prescribed from 
er by the simple rule (13.1011 . Each summand of (14.311 means 1 or 0 depending on whether such a 
boundary condition is satisfied or not for the chosen b 3 ® b 2 ® b 3 . 

Now we are in the position to convert this formula into a matrix product form by invoking 
Corollarv l2.8l To each vertex of (14.31) . substitute (12.211) or equivalently the graphical form (12.2211 of 
the combinatorial R. It amounts to putting the three “skewers” standing upright at the vertices. 
The result takes the form 


F(<T) = Tr(X ai ---X aL ), (4.4) 

where the trace extends over F® 3 and the factor X ak corresponds to the fcth layer of the skewers. It 
has the same structure as (14.31) but now each vertex represents an Ao-valued Boltzmann weight of 
the five-vertex model (12.2011 . Accordingly each arrow carries 0 or 1, and the sum over b 3 Cg) b 2 CS> b 3 
should be replaced by the configuration sum of the five-vertex model. What about the boundary 
condition on them? It should be imposed so as to reflect the table of Example 13.21 in the diagram 
(14.31) . In this way we find 




t 

o 

II 

M 

J l 


0 


0 





t 

II 

M 

J I 


0 


i 





t 

*2=£ 


1 


1 





t 

* 3 = £ 


i 


i 



(4.5) 

The sums here extend over all the configurations of the five-vertex model under the specified 
boundary conditions. The resulting objects take values in A® 3 . Explicitly they read 




+ 


+ 


+ 


+ 


t 

A 0 = 


+ 


+ 


+ 


+ 










= 1 <g> 1 ® 1 + a+ ® 1 ® 1 + k ® a + C3> 1 + a~ (g> a + ® a + + 1 <g> a + (g> a + 




+ 


+ 


t 

X, = 


+ 


+ 








= k<g>k<S> 1 +a <g)kig>a + +lcg}k(g>a + 




+ 


+ 


t 

a 2 = 


+ 


+ 








= l®a (gik + a+igia <g) k + k<g>l(g)k 
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+ 


t 


+ 


t 

*3 = 


+ 


+ 


+ 









= 1 ® a ® a + a + <g> a ® a + k ® 1 <g> a +a ® 1 ® 1 + 1 <g) 1 <g) 1 

Here the edges are colored black or red depending on whether they assume 0 or 1, and the three 
components in the tensor product are arranged so as to correspond to the vertices at (top left) ® 
(bottom left) ® (top right). See (12.2011 . 

Again general n case is similar and clear. By imposing the boundary condition according to 
(13.131) . the operator A, (0 < * < n) is given graphically as follows: 



(4.6) 


There are n incoming arrows from the left and n-outgoing ones from the top. It represents an 
element of 11 2 . The formula (j4.4|) remains the same if the trace is understood as the 

one over p®n(n i)/ 2 _ j n g enera i ft ft normalized so that P(cti, ..., <jl) = 1 if cti > ••• > cjl 
(which means ay = n, <jl = 0 as we are in a basic sector). As another example, the configurations 
contributing to X\ with n = 4 are as follows: 


^ ^ ^ ^ ^ ^ ^ ^ 


Accordingly we have X\ = kT-l-k-l-k+a~-a + -l-k-l-k+l-a + -l-k-l-k+l-k-Ta~-a + -k+a + -k-l-a~-a + -k+ 
kkl-la + k+la _ -a + -a _ a + -k+a + a^-a + a”a + k+ka~a + -la + k+a _ -l-a + -l-a + -k+l-l-a + -la + -k, 
where • denotes ® and A® 6 is ordered as (top line) ® (middle line) ® (bottom line) where within 
each line they are ordered from the left to the right. 

To summarize, our Xi is a corner transfer matrix of the q = 0-oscillator valued five-vertex model. 

The steady state probability P(er) (14.41) is a partition function of the 3D lattice model associated 
with the 3D L operator at q = 0 (12.201) . It is a 3D system having the prism shape region (14.61) xZ^ 
with the boundary condition on one of the surfaces specified by cr. The Xi plays the role of a 
layer-to-layer transfer matrix. In general it does not obey the recursion relation with respect to 
n of the form X-™' 1 = J2o<j<n-i ® ^ f° r some a ^ S Ajf" -1 . As a result it is different 

from those in mum except one particular case of n = 3 in the latter and the simplest nontrivial 
case n = 2: 

A 0 = l + a+, X\ = k, A 2 = l + a". 


5. Discussion 

5.1. Summary. In this paper we have revealed that the Ferrari-Martin algorithm for constructing 
the steady state of the n -TASEP in. is most naturally formulated in terms of a combinatorial R in 
crystal base theory. Combined with the factorized form of the combinatorial R originating in the 
tetrahedron equation (2Uj, it has lead to a new matrix product formula (14.41) for the steady state 
probability. Our operator Xi (14.61) is a corner transfer matrix of the g-oscillator valued five-vertex 
model at q = 0, which is a configuration sum in a cross-channel of the original problem. Whether 
such a result admits generalizations to similar systems like ASEP [22], open boundary conditions 
and the large list of factorized R matrices for other quantum groups [20] is a natural question to 
be investigated. 
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5.2. Generalized dynamics. Let us comment on a possible generalization of the time evolutions 
Ti and ts in Pronosition l3.4l from the viewpoint of crystal base theory. Let B l ’ s (1 < l < L, s > 1) be 
the crystal of the irreducible U q (sIl) module (so-called Kirillov-Reshetikhin module) corresponding 
to the l x s rectangular shape Young diagram mna. Elements of B l,s are labeled with the semi¬ 
standard tableaux on the Young diagram with entries {1,..., L}. The family B l,s includes B l (12.11) 
as B l}1 = B l , and again there is a bijection B l,s g) B l ,s —>• B l ,s g) B l ' s called the combinatorial 
R. For an arbitrary pair of (l, s ) and (V : s'), an algorithm for finding the image of the R is known 
[23] which reduces to the NY-rule for s = s' = 1. 

Associated with any crystal B r,s we introduce the maps T = T r,s : B{m) ® B r ’ s — > B r ’ s g>£>(m) 
and T u '■ £>(m) —»■ B{ m) for each u G B r ’ s by a scheme similar to (13.141) and (13.161) . Namely we set 
T(bx ®— g) b„ ® u) = u' g) g) • • • g) b' n and T U {B\ g> - - - <8> b n ) = b^ g> • • • g) b(,, where u! G B r ' s 
and b' G B li are obtained by sending u G B r ’ s to the left through tq g) • • • g) b„ by a successive 
application of the combinatorial R: B li g) B r ’ s —>• B r,s g) B li . Maps like T u have been studied 
extensively as a time evolution with carrier u in the context of integrable cellular automata known 
as box-ball systems and their generalizations M- Based on computer experiments we propose 

Conjecture 5.1. For cr G S(m), let n~ 1 (cr) C B{ m) denote the preimage of the map n II 3.2(A) . 
Then the set tt o 7) l (7r _1 (<x)) consists of a single element for any u G B r,s . 

Admitting the conjecture we are able to introduce a new time evolution t u on S'(m) by t u := 
7r o Tu o 7T ~ 1 . The definition is equivalent to replacing Ti and r,; in the commutative diagram of 
Proposition l3.4l bv T u and t u , respectively. The detail of the dynamics will be illustrated elsewhere. 
The map T = T r,s is invertible as it is a composition of combinatorial R's. 

Now we introduce new stochastic processes on B( m) and 5(m) analogous to the multiline process 
and the n-TASEP, respectively. They are defined by letting each state undergo the respective time 
evolutions T u and t u with an equal probability with respect to u G B r,s . Then from the invertibility 
of T, exactly the same argument as [TT] p821] proves that the steady states are given by the uniform 
probability distribution for the former and by (13.211) for the latter. In this way we find a curious 
feature of the probability distribution (13.211) that it yields a simultaneous steady state for the family 
of stochastic processes on 5(m) associated with the crystals B r,s as well as for the n-TASEP. 


5.3. Hat relations. Our derivation of (14.41) and (14.61) in this paper is based on Proposition 
13.51 which is attributed to m However it is also possible to establish them directly follow¬ 
ing the idea in [Tj. To explain it let us write the action of the local Hamiltonian (13.41) as 
h\a,P) = Eo< 7 ,5<n^l7^)' Then for the quantity Tr(A CTl ■ ■ ■ X aL ) in (|ili|) to give the steady 
state probability, it is sufficient to construct the operators X 0 ,..., X n G End(F®"( n-1 )/ 2 ) satisfy¬ 
ing the so called hat relation [7]: 

X a Xp — X a Xp = (hXX) a< p := Y, Ks x i x 5 (0 <a,p<n). (5.1) 

0<7,<5<n 

Theorem 5.2. Up to a freedom Xf —> Xi + const Xi at least, such an Xj is realized as 

OL i 0)2 Oin 


Xi — + • • • + 0! n ) 


(5.2) 

where the sum ranges over the same set of configurations of the five-vertex model as in Xi B- 
The only difference is the extra coefficient a\ + ■ ■ ■ + a n depending on the n outgoing arrows 

cii = 0,1. 
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For instance one has Xo = a + 0 1 (g) 1 + k 0 a + (^) 1 + a - (g) a + a + + 2(1 0 a + 0 a + ) for 
n = 3 from the graphical representation of Xo in Section 14.21 The hat relation ()5.1| ) turns out to 
be a consequence of its far-reaching generalization into a 3D lattice model obeying the tetrahedron 
equation. An exposition of the detail is beyond the scope of this paper and will be presented in 
[T9] together with a proof of Theorem 15.21 

Acknowledgments 

The authors thank Yasuhiko Yamada for useful discussion and Chikashi Arita for a communi¬ 
cation. A. K. thanks Tetsuo Deguchi, Akinori Nishino, Tomohiro Sasamoto and Hal Tasaki for 
kind interest and valuable comments at Mathematical and Statistical Physics Seminar at Ochano- 
mizu University, 6 June 2015, where a part of the work was presented. This work is supported by 
Grants-in-Aid for Scientific Research No. 15K04892, No. 15K13429 and No. 23340007 from JSPS. 

References 

[1] C. Arita, A. Ayyer, K. Mallick, S. Prolhac, Recursive structures in the multispecies TASEP J. Phys. A: Math. 
Theor. 44 (2011) 335004 (26pp). 

[2] C. Arita, A. Kuniba, K. Sakai and T. Sawabe, Spectrum in multi-species asymmetric simple exclusion process 
on a ring. J. Phys. A: Math. Theor. 42 (2009) 345002 (41pp). 

[3] R. J. Baxter, Exactly solved models in statistical mechanics , Dover (2007). 

[4] V. V. Bazhanov and S. M. Sergeev, Zamolodchikov’s tetrahedron equation and hidden structure of quantum 
groups, J. Phys. A: Math. Theor. 39 3295-3310 (2006). 

[5] R. A. Blythe and M. R. Evans, Nonequilibrium steady states of matrix product form: A solver’s guide, J. Phys. 
A: Math. Theor. 40, R333 (2007). 

[6] E. Date and M. Okado, Calculation of excitation spectra of the spin model related with the vector representation 
of the quantized affine algebra of type Int. J. Mod. Phys. A 09 399-417 (1994). 

[7] B. Derrida, M. R. Evans, V. Hakim and V. Pasquier, Exact solution of a ID asymmetric exclusion model using 
a matrix formulation, J. Phys. A: Math.Gen. 26 1493-1517 (1993). 

[8] V. G. Drinfeld, Quantum groups , In Proceedings of the International Congress of Mathematicians, Vol. 1, 2 
(Berkeley, Calif., 1986), pp798-820, Amer. Math. Soc., Providence, RI, (1987). 

[9] V. G. Drinfeld, On some unsolved problems in quantum group theory. In “Quantum groups”, Lect. Notes in 
Math. 1510 (1992) pl-8. 

[10] M. R. Evans, P. A. Ferrari and K. Mallick, Matrix representation of the stationary measure for the multispecies 
TASEP, J. Stat. Phys. 135 217-239 (2009). 

[11] P. A. Ferrari and J. B. Martin, Stationary distributions of multi-type totally asymmetric exclusion processes, 
Ann. Probab. 35 (2007) 807-832. 

[12] G. Hatayama, A. Kuniba, M. Okado, T. Takagi and Z. Tsuboi, Paths, crystals and Fermionic formulae, Prog, 
in Math. Phys. 23 Birkhauser (2002) p205-p272. 

[13] J. Hong and S-J. Kang, Introduction to quantum groups and crystal bases, Graduate Studies in Math. 42 
(AMS, 2002). 

[14] R. Inoue, A. Kuniba and T. Takagi, Integrable structure of box-ball systems: crystal, Bethe ansatz, ultradis¬ 
cretization and tropical geometry, J. Phys. A. Math. Theor. 45 (2012) 073001 (64pp). 

[15] M. Jimbo, A q- difference analogue of U(q) and the Yang-Baxter equation, Lett. Math. Phys. 10 (1985) 63-69. 

[16] M. Kashiwara, On crystal bases of ^-analogue of universal enveloping algebras, Duke Math. J. 63 (1991) 
465-516. 

[17] S-J. Kang, M. Kashiwara, K. C. Misra, T. Miwa, T. Nakashima and A. Nakayashiki, Affine crystals and vertex 
models, Int. J. Mod. Phys. A 7 (suppl. 1A), (1992) 449-484. 

[18] S-J. Kang, M. Kashiwara, K. C. Misra, T. Miwa, T. Nakashima and A. Nakayashiki, Perfect crystals of quantum 
affine Lie algebras, Duke Math. J. 68 (1992) 499-607. 

[19] A. Kuniba, S. Maruyama and M. Okado, Multispecies TASEP and the tetrahedron equation, in preparation. 

[20] A. Kuniba, M. Okado and S. Sergeev, Tetrahedron equation and generalized quantum groups, J. Phys. A: 
Math. Theor. 48 (2015) 304001 (38pp). 

[21] A. Nakayashiki and Y. Yamada, Kostka polynomials and energy functions in solvable lattice models, Selecta 
Mathematica, New Ser. 3 (1997) 547-599. 

[22] S. Prolhac, M. R. Evans and K. Mallick, The matrix product solution of the multispecies partially asymmetric 
exclusion process, J. Phys. A: Math. Theor. 42 (2009) 165004 (25pp). 

[23] M. Shimozono, Affine type A crystal structure on tensor products of rectangles, Demazure characters, and 
nilpotent varieties, J. Alg. Comb. 15 (2002) 151-187. 

[24] A. B. Zamolodchikov, Tetrahedra equations and integrable systems in three-dimensional space, Soviet Phys. 
JETP 79 641-664 (1980). 


16 


ATSUO KUNIBA, SHOUYA MARUYAMA, AND MASATO OKADO 


E-mail address-. atsuoSgokutan. c .u-tokyo.ac. jp 

Institute of Physics, University of Tokyo, Komaba, Tokyo 153 - 8902 , Japan 
E-mail address: maruyamaSgokutan.c.u-tokyo.ac.jp 

Institute of Physics, University of Tokyo, Komaba, Tokyo 153 - 8902 , Japan 
E-mail address: okadoSsci.osaka-cu.ac.jp 

Department of Mathematics, Osaka City University, 3 - 3 - 138 , Sugimoto, Sumiyoshi-ku, Osaka, 558 - 
8585 , Japan 


